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Abstract 

We investigate self-similar solutions of evolution equation of a (1+1)- 
dimensional field model with the V-shaped potential U(ip) = \<p\, where ip 
is a real scalar field. The equation contains a nonlinear term of the form 
sign(ip), and it possesses a scaling symmetry. It turns out that there are 
several families of the self-similar solutions with qualitatively different be- 
haviour. We also discuss a rather interesting example of evolution with non 
self-similar initial data - the corresponding solution contains a self-similar 
component. 
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1 Introduction 




sign((p(x,t)) 



(1) 



where <p is a real scalar field, and x, t are dimensionless variables obtained by 
appropriate rescalings of the physical position and time coordinates. It follows 
from the Lagrangian 



where // = 0, 1, d = d t , d\ = d x . The sign function has the values ±1 when 
(p ^ 0, and if <p = 0. Because of the sign(tp) term, Eq. (1) is nonlinear in 
a rather interesting way. The corresponding field potential U(<p) = \<p\ has the 
minimum at tp = 0, and the field tp can oscillate around the equilibrium value 
(f = 0, but Eq. (1) cannot be linearized even if amplitude of the oscillations is 
arbitrarily small. 

The signum-Gordon model has a rather sound physical justification. It has 
been obtained in a continuum limit (or in a large wavelength approximation) to 
certain easy to built mechanical system with a large number of degrees of freedom, 
see fTJ [2) for details. The model can also be applied in a description of pinning 
of an elastic string, which represents a vortex, to a rectilinear impurity in the case 
both the string and impurity lie in one plane. Yet another interesting application 
is to dynamics of a system of two global strings on a plane. Such strings are 
represented by the real functions ipi(x, t), ip2(x, t), and in certain approximations 
their dynamics is summarized in the following Lagrangian 



where V is an extrapolation of the well-known logarithmic interaction potential 
between separated global strings to small values of \ipi — ip 2 \- 



where a is a positive constant (see 01 for a detailed discussion of the inter-string 
potential). For small |-0i — i/^l/ we nave 

V(V>i -^2) « |^i -fal- 
'We thank Benny Lautrup from NORDITA for suggesting this name. 



L = -d^(p&*(p - \<p 
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In this case the Euler-Lagrange equations for ipi, ijj 2 obtained from Lagrangian L s 
imply that ip = \{$\ — ^2) obeys the signum-Gordon equation. The solutions 
presented in our paper describe particular cases of time evolution of all those 
systems. The present paper, however, is focused on self-similar solutions of the 
signum-Gordon equation rather than on its applications. 

The signum-Gordon equation possesses the exact scaling symmetry: if (p(x, t) 
is its solution then 

<px(x,t) = \ 2 <p(j,j) (2) 

is a solution too, for any constant A > H This symmetry is one of the most 
interesting features of the model. As always when there is a symmetry, one may 
search for solutions which are invariant under the symmetry transformations. In 
the case of scaling they are the self-similar ones. In paper |Z| an example of such 
solutions of the signum-Gordon equation has been given. Self-similar solutions of 
nonlinear equations play important role in nonlinear dynamics. They have plenty 
applications - it is impossible to list them here. Beautiful presentation of the self- 
similarity with its applications is given in [4J. A shorter introdunction to this topic 
can be found in 0. 

The present paper is devoted to a thorough investigation of self-similar solu- 
tions of the signum-Gordon equation. Our main finding is that, rather surpris- 
ingly, there are quite many classes of such solutions. The solutions presented 
in B21 form merely a measure zero subset in just one such class: the segment 
So = 0, —1/2 < So < 1/2 in Fig. 1. The various types of solutions differ in 
particular by the number of isolated zeros of the field cp: they can have infinitely 
many, just one or no isolated zeros. Funnily enough, all self-similar solutions are 
composed of quadratic polynomials in x and t. 

We have also investigated the evolution from certain particular initial data 
which are not self-similar, with d x cp discontinuous at certain point. While the 
late time behaviour of the corresponding solution remains mysterious, we have 
noticed a very interesting phenomenon which occurs during the early stages of 
the evolution. At first, the point of the discontinuity of d x ip moves with a con- 
stant velocity equal to —1 ('the velocity of light') until d x ip> becomes a continuous 
function of x. From this moment on, the evolution of cp locally, in a vicinity of 
the former point of discontinuity, follows a self-similar solution with one isolated 
zero. This behaviour is universal in the sense that it does not depend on details of 

2 Transformations of the form (2) with A < are obtained as a product of the scaling trans- 
formation (2) with the reflections x — > —x, t — > —t. Eq. (1) is invariant with respect to such 
reflections, and also with respect to 1+1 dimensional Poincare transformations. 
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the initial data. Solutions of this class might be relevant for the description of the 
process of coalescence of two strings, see Fig. 9. 

All our analytic solutions have been checked by comparisons with purely nu- 
merical solutions of Eq. (1). Actually, there was a very fruitful interplay of the 
two methods: certain solutions were seen first numerically and later obtained in 
the analytic form, whereas with others it was the other way round. 

The plan of our paper is as follows. In Section 2 we discuss the Ansatz and 
initial data for the self- similar solutions, and we present a map of such solutions. 
Section 3 is devoted to a detailed presentation of all classes of the self-similar 
solutions: in subsections 3.1-3.5 we give explicit forms of all solutions. In Section 
4 we consider the initial stages of evolution in the particular case of non self- 
similar initial data. Section 5 contains a summary and remarks. 

2 The Ansatz and initial data 

The physical context in which the signum-Gordon equation (1) appears implies 
that the most interesting solutions are such that ip(x, t) and their first order partial 
derivatives d t ip, d x tp are continuous functions on the (x, t) plane, perhaps except 
across the lines X — Xq ± t (the characteristics). For example, a discontinuity of 
d t tp as a function of time t at a certain instant ti would mean that the velocity 
of a certain part of a continuous mechanical system (e.g., of a string) suddenly 
jumps in spite of the fact that the force, given by the sign(Lp) term, is always finite. 
Obviously, such discontinuities are unphysical. On the other hand, Eq. (1) implies 
that at least one of the second order derivatives d^(p, d%(f has to be discontinuous 
when the piecewise constant function F(x,t) = sign(ip(x,t)) changes its value 
from to ±1. At such a point ip = 0, and in a certain neighbourhood of it if > 
or (p < 0. Then, the left and right second derivatives of (p can have different 
values - if this is the case then the ordinary second derivative of if does not exist, 
strictly speaking. Nevertheless, such solutions are physically admissible. The 
proper mathematical framework for discussing them is based on the notion of 
weak solutions, see, e.g. J6J. 

In the case of transformation (2) the scale invariant Ansatz can be taken in the 
form 

ip(x,t) = x 2 S(y), y=-. (3) 

x 

We will consider solutions of Eq. (1) for t > with initial data specified at t = 0. 
Solutions for t < can be easily obtained with the help of the time reflection 
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symmetry. It will turn out that S(y) can be a quadratic, linear or constant function 
of y. Therefore, ip will be continuous at x = in spite of the fact that the scale 
invariant variable y is singular at that point. The equivalent Ansatz 

V (x,t)=t 2 T(^) 

avoids the superficial singularity at x = 0, but it is less convenient for incorporat- 
ing the initial data. 

The self-similar initial data have the form 

, , f R x 2 for x < 0, , / R x for x < 0, 

^°) = { S x* for *>0, MM)U = | s x for x > 0, (4) 

where i? , -Ro, 5*0, 5*0 are constants. Such initial data are self-similar because the 
point x = is not shifted by the rescaling x — > x/A. In general, such </?(x, 0) is of 
the C^-R) class, while d t (p(x,t)\ t=Q is of the class C°(R). 

We shall restrict our considerations to the case R = R = 0. The other 
choice: S — S — and R , R arbitrary is related to the previous one by the 
spatial reflection x — > —x, which is a symmetry of Eq. (1). Thus, in the main part 
of our paper we assume that 

( n , j for x < 0, a , f for or < 0, 

^ (X ' 0) = \ 5ox 2 for x > 0, M^OU = { ^ for x > 0. (5) 

Such initial data correspond to the following conditions for the function S(y): 

hm + S(y) = S , \\m_S(y) = 0, ^lim S'(y) = S , \im_ S'(y) = 0. (6) 

Let us remark that solutions for more general self- similar initial data (4) can be 
obtained in certain cases just by combining solutions obeying (5) with the ones 
obtained from them by applying the spatial reflection, but there are also cases in 
which such approach does not work. 

The Ansatz (3) reduces Eq. (1) to the following ordinary differential equation 
for the function S(y) 

(1 _ ^ S " + 2yS'-2S = -sign(S), (7) 

where S' = dS/dy, and the sign(S) function has the values +1 or —1 for S > 
or S < 0, respectively, and sign(0) = 0. Notice that at the points y — ±1 
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(which correspond to characteristics of Eq. (1)) the coefficient in front of the 
second derivative term in Eq. (7) vanishes. This has the consequence that the first 
derivative dS/dy of the solution does not have to be continuous at these points. 
Below we shall make use of this possibility. In order to obtain weak solutions of 
this equation we first solve it assuming that S > or S < or S = 0. Such partial 
solutions have a rather simple form of quadratic polynomials in y: 

when S>0: S+(y) = -^{y 2 + 1) + \y + \ 

when S<0: S_(y) = |(y 2 + 1) - jy - \ 

where a, (3, a', (3' are arbitrary constants. The polynomial solutions are valid on 
appropriate intervals of the y-axis determined by the conditions S > or S < 0, 
respectively. There is also the trivial solution 

S (y) = o. 

Next, we match such partial solutions requiring that S and S' are continuous 
functions of y, except at the point y — 1, where only continuity of S is required. 
It turns out that the continuity of only S at y — 1 is sufficient to ensure the implied 
by wave equation (1) continuity of 

d(p 1 (dip d(p\ 
dx + 2 \dx dt J ' 

where + t, across the characteristic line x = t. The matching condi- 

tions together with the initial data (5) determine the constants a, (3, a', (3' from 
the partial solutions. In this manner we have constructed exact, self-similar weak 
solutions of the Cauchy problem (5) for all values of parameters So, So- Because 
x E (— oo, oo), the variable y can take all real values. The factor x 2 in (3) makes 
the solutions finite at x = 0. 

Depending on the values of So, So the self-similar solutions can have quite 
different forms. As one can see from the 'map' presented in Fig. 1, the space of 
such solutions of the signum-Gordon equation is surprisingly rich. 
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Figure 1: The map of the self- similar solutions of the signum-Gordon equation 
(1). Each point in the (So, So) plane represents one solution. The symmetry of 
the picture with respect to the reflection (S , S ) — ► (—So, — So) is related to the 
reflection symmetry cp — > —if of the signum-Gordon equation. Because of this 
symmetry it is sufficient to discuss solutions from the half-plane So > 0. The 
dashed lines are auxilliary: two coordinate lines, and the two lines S = — 2S 
which are asymptotically tangent to the lines v — 0. The continuous lines are 
border lines between classes of solutions, but their points also represent some 
special solutions discussed in subsections 3.1 (the v = 0, So > rectilinear half- 
line), 3.2 (the v = 0, So > curved line), 3.4 (the open segment connecting the 
point (0,-1/2) with the point (1/4, 0)), and 3.5 (the v = 1, S > half-line), vis 
a certain velocity defined in the text. The lower, curved v = line is a cubic curve 
given by equation (16), the other one is obtained by the reflection. The rounded 
brackets, which enclose two line segments - the edges of the rhomb III, indicate 
that these segments do not contain ends - the end points belong to the two v — 1 
half-lines. 
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3 Self-similar solutions 



Below we give analytic forms of the self-similar solutions with initial data (5). We 
assume that S > 0. Solutions with negative values of S can be obtained by the 
reflection <f)(x, t) — > —(/>(x, t), which in particular means that S — > — So, So - ► 



3.1 Solutions of types: la, Ha, v = 

These solutions lie in the open region above the rhomb III and above the v = 
1, So > line, see Fig. 1. The solutions are obtained by matching the trivial 
solution with the partial solution S + at the point y — 1/v which corresponds to 
the line x = vt. The velocity v has values in the interval (—1,1), and it is given by 
formula (8). It turns out that one has to use a second solution of the type S+ which 
matches the previous one at the point y — 1 corresponding to the characteristic 
line x = t. That latter S + solution has the form S + S y — (1/2 — S )y 2 , where 
So > 0, S are the parameters specifying the initial data according to formula (5). 
The velocity v is given by the formula 



which follows from the matching condition at the point y — 1. In the region la the 
velocity has values from the interval (—1, 0), while in the region Ha v € (0, 1). On 
the line 25*0 + So — 1, which separates the two regions, v — 0. The full solution 
has the form 



-S . 




for 
for 

\)t 2 for 



x < vt, 

vt < X < t, 
x>t. 



(9) 



It is depicted in Fig. 2. 



The half-line 2S + S = 1, S >0,v 
Ha, includes the static solution 



0, which separates the classes la and 




for x < 0, 
for x > 0, 



(10) 
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A 



<Po t 

x = vt 

Figure 2: Picture of a solution of type la for which — 1 < t> < 0. f is the velocity 
of the point on the x axis at which cpi merges with the trivial solution tp = 0. The 
two vertical arrows indicate that values of ip increase with time - the cp > part 
moves to the left. Solutions of the type Ha have a similar shape, but they move to 
the right (0 < v < 1). 

for which So = 1/2, So = 0. Other solutions from that line are time dependent. 
They coincide with the static one in the growing with time interval < x < t, 
while at the points x > t they have the form ip 2 {x, t), see the second and third 
lines in formula (9). Asymptotically, as t grows to infinity, such cp entirely covers 
the static solution. 

The signum-Gordon equation is invariant under Lorentz transformations 

i! = j(t — wx), x' = j(x — wt), 

where 7 = 1/yT — w 2 , \w\ < 1. Such transformations preserve self-similarity of 
solutions. One can see that the static solution (10) is transformed into the solution 
with So = 7 2 /2, So = —wy 2 . In this case formula (8) gives v = w. For these 
particular solutions the functions cpi and (p 2 in formula (9) have the same form. 

3.2 Solutions of types: lb, lib, v = 0. 

These solutions lie in the open region below the rhomb III and the v = 1, S > 
line, see Fig. 1. Solutions of this type have the form shown in Fig. 3. Solutions 
from the half-line 5*0 = 0, So < —1/2 are not considered here because they can 




9 



be obtained from solutions lying on the S = 0, S > 1/2 half-line with the help 
of the reflection 0^—0. 




Figure 3: Picture of a solution of type lb for which — 1 < v < 0. The arrows 
indicate that values of decrease with time - the ip < part expands in both 
directions. Solutions of the type lib have a similar shape, but they move to the 
right (0 < v < 1). 



They are composed of the trivial solution S = in the region y > 1/v > 1, 
one solution S + in the region y < l/v\ < 1, and two solutions S- in the region 
1/vx < y < 1/v which continuously match each other at the point y = 1. The 
pertinent formula for ip has the following form 



(p(x,t) 



<Pg 
<Pi 

^3 





(x-vt) 2 
2(l-v 2 ) 

_|_ f2\ _ a2^,f _ 1^,2 

5 x 2 + S tx + (S - \)t 2 



for 
for 
for 
for 



x < vt, 

vt < X <t, 
t < X < Vit, 
X > Vit. 



(11) 



Here S , S are given in the initial data (5). The velocity v% is obtained from the 
condition ip 3 (vit, t) = 0, and a 2 , (3 2 from the matching conditions at x = V\t, i.e. 
at y = 1/vi < 1 . Straightforward calculations give 



— So 



C2 

~2S 



ASl + 2So 



(12) 
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Vl + (l +v ^Si-ASi + 2S 
«2 = 2 1 — , (13) 



ft = \. 2 ; " -. (14) 



v 2 1 +2v 1 ^ /Si-A Si + 2So 
v\-l 

Finally, the velocity v is obtained from the condition of continuity of ip{x, t) at the 
point x — t : 

2j3 2 - a 2 



(15) 



" 2 + a 2 - 2[3 2 ' 

Formulas (12-15) imply that t> x > 1 and — 1 < v < 1. 
The condition t> = is satisfied by (So, S ) such that 

' = \ + w <16) 

where 

^ = 250-50, C = 2S + S . 

These points form the lower curved line in Fig. 1. Solutions of this kind are 
related to another static solution of Eq. (1), namely to 

, , f for x < 0, 

= { -x*/2 for x > 0, (17) 

for which S = —1/2, S = 0. These v — solutions are time dependent. 
Their </?i part, see formula (11), coincides with the static one in the growing 
with time interval < x < t, while at the points x > t their form is given by 
(p 2 (x, t), y?3(x, t). Asymptotically, as t grows to the infinity, we recover the static 
solution. 



3.3 Solutions of type III 

Let us try to put together a number of the partial solutions S + , S_ in order to 
cover as large as possible interval of the y > half-axis, see Fig. 4, where 
Sk, k — 1,2,..., denote the consecutive partial solutions S±. 

The restriction to y > 0, see Fig. 4, is quite natural. The point is that y = 
0_ corresponds to x = — oo, while x = + to x = +oo. Hence, there is no 
physical reason for demanding that S(y) is continuous at y = - the continuity at 
y — would mean that we arbitrarily impose the restriction y?(oo, t) = ip(— oo, t). 
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Therefore, the solution Si is terminated at y = 0, and in the region y < we take 
the trivial solution S (y) = 0. On the other hand, both y = — oo and y = +oo 
correspond to x — 0, where the solution cp is required to be a continuous and 
differentiable function of x. 




Let us write the partial solutions in the form 

S k (y)= 1 -(-l) k [(3 k (y 2 + l)-a k y-l}. 

Equation (1) implies that the following matching conditions at the points y = 
Ofe, k > 1, have to be satisfied: 

S k {a k ) = = S k+ i(a k ), S' k (a k ) = S' k+1 (a k ), (18) 

provided that a k ^ ±1. When y = ±1, only continuity of S(y) has to be required. 
The matching conditions (18) yield the recurrence relations: 

4a fc 2 _ 2a k - (1 + a|)a fe _i 

1 — 1 — a k 1 + a k — 2a k a k -i 

(19) 

where aii = ai > 0, a < 0, and 

/3i = , a x = . (20) 

1 — a^ai 1 — a ai 
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The relations (20) follow from the fact that a , a\ are the zeros of the polynomial 
S\(y). It is clear from Fig. 4 that a < 0, a\ > 0. The polynomial S\(y) has 
zeros at y = ao,a±. 

The solution is single- valued when a k +i > a k . Simple calculations show that 
this condition taken for k — 1, 2 implies that ai < 1, a 2 < 1, and that a > — 1. 
Furthermore, we notice that a = —1 or ai = 1, always give a 2 = 1, what implies 
that a k = 1 for all k > 2. Such solutions, i.e., with a 2 = 1, belong to the classes 
of solutions discussed in subsections 3.4, 3.5, except for the case a = — 1, ai = 
in which we obtain the trivial solution Sq = 0. Therefore, in the remaining part of 
this subsection we assume that 



-1 < a < 0, < ai < 1. 



Recurrence relations (19) have the following general solutions: 



1 



1 + r 



1 



< p r f 



Y — pr k + (— 1' 



v 



-i)V 



(21) 



and 



where 



1 

2 + 



2(1 



pr 



k-l 



k , 

pr + 



a k 



1 — ai 



P 



-l) k (l 



1 - q k 



p 



q) 



pk 1 _|_ qk 



■l) k Pi, (22) 



(23) 



P 



a 



1 + a 



Q 
P 



(24) 



1 + Oi ' 

Notice that < q < 1 and p > 1. 

Solutions (21-23) have been conjectured after seeing several first iterations 
of recurrence relations (19), and checked by substitution into these relations. In 
paper j6l we were able to find only solutions with a = 0, a\ < 1. They lie on 
the segment 5 = 0, < S < 1/2 inside the rhomb III in Fig. 1. 

Formula (23) implies that a k < 1, and that — > 1 when k — >• oo. There- 
fore, the piecewise polynomial solution constructed above covers only the inter- 
val [0, 1). For this reason we introduce a special notation for it, namely S pp (y). 
Calculations show that S pp (y) vanishes when y — > 1. On the other hand, the first 
derivative S' Ay) does not vanish in that limit, but it remains finite. 

This solution can be extended to the full range of y. We have already taken the 
trivial solution So(y) = in the region y < 0. In the region y > 1 we also take 
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the trivial solution which is consistent with continuity of S(y) at y = 1. The first 
derivative dS/dy in general is not continuous at y = 1, but this is not forbidden by 
Eq. (7) because the coefficient 1 — y 2 in front of 5" vanishes when y — 1. Notice 
that y = ±1 correspond to the characteristics a; = ±t of the signum-Gordon 
equation. 

The solution we have just constructed has the following form: 

f A - / = for x <t, 

_ \ Vpp = x 2 S pp {{) for x>t. 

Snapshot of the solution is presented in Fig. 5. 




Figure 5: The self-similar solution of type III. The zeros x fc = v^t of (p accu- 
mulate at x = t because Vk = l/cife, where a k — > 1 when k — > oo. The arrow 
indicates that whole structure moves to the right. 



The zeros of are given by the formula Xk(t) = t/a^. They move with the 
constant velocities v^. = l/a^ > 1 . 

The parameters So, So, which specify initial data (5) for these solutions, enter 
the first polynomial Si . Hence, 

S = Si(0) = i(l - ft), So = S[(0) = ^, (26) 

where ai, Pi have the form (20). Because a , a x are zeros of this polynomial, 
they are related to the initial data by the following formulas 

So - Va So + VA 



1 - 2S 1-2S 
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where A = Sq — AS^ + 2S . By the assumption, for the considered solutions 
-1 < a < 0, < ai < 1, see Fig. 4. It follows that 

2S -^<S <^-2S , 0<S <^. (27) 

These inequalities correspond to the So > part of the interior of rhomb III in 
Fig. 1. 



3.4 Solutions from the bottom-right open edge of rhomb III 

In the present case we take the trivial solution in the regions y < and y > 1 . At 
the point y — 1 the trivial solution matches the partial solution S_. Because that 
point lies on the characteristic line of the signum-Gordon equation, we demand 
only continuity of S(y). The solution S_ in turn matches a partial solution S + 
at a point y x from the interval < y\ < 1. Here we also demand continuity 
of derivatives dS/dy. Of course, both S-(yi) and S+(yi) are equal to 0. The 
calculations are straighforward. It turns out that 

= 4 " L <28) 



This formula implies that v x > 1. The solution (p has the following form 

<p(x,t) = < 



(p — for x < t, 

¥1 = 2 (vi-i) ( x -t)( x - v it) for t<X<Vit, (29) 



</? 2 = (x + ^[^(x + t) — |] for x>vit, 



see Fig. 6. 



These solutions have the following initial data 

S = 2S -^, 0<S <^. (30) 

They lie on the open segment which is the interior of the > 0, S < edge of 
rhomb III in Fig. 1. 

Let us note that these solutions can be regarded as a limiting case of solutions 
discussed in the previous subsection obtained by putting a = — 1 and keeping a\ 
in the open interval (0, 1). 



15 




Figure 6: The self- similar solution (29). 

3.5 Solutions with v = 1 

In this case we combine the trivial solution and an S + solution. They match 
continuously at the point y — 1. 



<p(x,t) 



y?o = for x < t, 

<p 1 = (x-t)[S (x-t)+t/2] for x>t, 



(3D 



see Fig. 7. In particular, for So = 1/4 we have ip(x, t) = (x 2 — t 2 )/A for x > t, 
and if = for x < t. 



A<p 



x = t 



Figure 7 : The self- similar solution (31). 
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The parameters So, S for these solutions obey the relation 

So = -2S + l -, (32) 

where S > 0. 

Some of these solutions are related to the ones discussed in subsection 3.3: if 
we take a\ = 1 and — 1 < a < 0, then only the first polynomial Si is present and 
it leads to the solution of the form (31). However, the restriction — 1 < a < 
implies that < S < 1/4 because S = — a /2(l — a ). This corresponds to that 
part of the v — 1, S > half-line, see Fig. 1, which is the upper-right edge of 
rhomb III. 

4 Non self-similar solutions with a self-similar com- 
ponent 

In this Section we would like to present a certain interesting class of solutions of 
the signum-Gordon equation which are not self- similar. The corresponding initial 
data for them are different from (5): they have the form 

**•") = { L t III W*.*)U=0, 03) 

where a > is a free parameter. The solution if was first obtained numerically 
in a context which is not relevant here. Computer simulations showed that, rather 
surprisingly, certain important quantitative characteristics of the evolution of ip 
apparently did not depend on the parameter a at all. Such numerical findings have 
motivated us to study these solutions in more detail. It has turned out that the 
solutions have a self-similar component during a certain finite time interval. For 
this reason we would like to present them here. 

Because the signum term in Eq. (1) remains constant until (p becomes equal to 
zero, the analytic solution of Eq. (1) with the initial data (33) can be constructed 
piecewise. In each interval on the x axis such that (p has a constant sign in it, 
one can use the well-known formula for general solution of the one-dimensional 
wave equation, suitably modified to incorporate the constant +1 (or —1) term in 
the equation: 

ip(x, i) = h(x — t) + g(x + t) + Cit 2 — c 2 x 2 , 
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where 2(ci + c 2 ) = ±1(= sign(ip)). The functions h, g and the constants Ci,c 2 
are determined from the initial conditions and from the matching conditions. The 
resulting solution has a relatively simple form when the time t is not too large. The 
evolution of ip can be divided into distinct stages. In the first one time t changes 
from to a. Then 



<p {x,t) = 

<P<a(x,t) 







ax — 



t)(4a + x-3t) 

+2 



for x < —t, 

for -t < x < t, (34) 

for x > t. 



The function ip a i remains positive in the interval x € (—t, t] until t — a. Notice 
that the border point betwen the supports of functions ip Q , ip al moves with the 
velocity —1 which does not depend on the slope a of the initial shape (33) of (p. 
This fact is related to the jump of the value of d x (p a at the point x = — t - such 
discontinuity is allowed for by the wave equation only on characteristics. Hence 
the border point between the supports of the functions (p , (p a i has to move with 
velocity -1 until d x (p al \ x= _ t vanishes. This happens at the moment t — a. 
At that moment 



(p a (x,a) 







a(x 



and 



(d t ip a )(x,a) 



-\{x + a) 
-a 



for x < —a, 
for —a < x < a, 
for x > a, 



for x < —a, 
for —a < x < a, 
for x > a. 



(35) 



(36) 



These values of </?, d t ip constitute the initial data for the second stage of evolu- 
tion of ip which lasts until t = 2a. Here the crucial observation is that (p a , d t ip a 
at the time a have the shape which coincides with the self- similar initial data 
(5) in which x is replaced by £ = x + a. The corresponding solution belongs 
to the class discussed in subsection 3.4. It has S = 1/8, S = —1/4, and 
v 1 = 3. It is interesting that the velocity v 1 does not depend on the parameter 
a. The y?i part, see Fig. 6, begins to appear at the point x = —a at the time 
t — a. It expands to the right with velocity +3, while the left end of its support 
moves slower, with velocity +1. Thus, we expect that in the interval x G (— oo, t] 
the solution has the form of time and space translated self-similar solution (29) 
(£ — > r = i — a, x £ = x + a). At the point x = t it continuously matches the 
function (p a2 from formula (34). With these hindsights one can easily construct 
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the analytic solution. In the time interval t E [a, 2a] we have (p(x, t) = (fb(x, t), 
where 



(p b (x,t) 



Here 





fbi = \{x - x (t))(x - 3f + 4a) 



Vb2 



[X 



t)(x - 3t + 4a) 



for 
for 
for 
for 



x < t — 2a, 
t-2a<x<3t-4a, 
3t - 4a < x < t, 
x > t. 

(37) 



xo(t) = t — 2a. 

Notice that (p b3 describes the "freely falling" halfline if = ax, which is the 
remnant of the initial data (33). 

At the time t — 2a another interesting thing happens. At this moment the sup- 
port of (fb2 becomes reduced to the point x = 2a, and simultaneously </?&3 (2a, 2a) = 
0. The right end of the support of (fbi, which moves with the velocity +3, hits the 
support of (fib3- At this time the length of the support of (fbi is equal to 2a. 

In the next time interval, t G [2a, 3a], the solution ip = <p c has the following 
form, see also Fig. 8, 



ip c (x,t) 



where 



' (fo = 

<Pci = ~\{x - x (t))(3t - x - 4a) 
V?e3 



1+2 



k y? c4 = ax - id 



for 
for 
for 
for 
for 



x < t — 2a, 

t -2a < x <4a 

4a - t < x < t, 



t, 



(38) 



3 2 3 7 o 2 1 7 \f0L m/ \ 

<£c2 = -x — ~.xt H — t +a at H — ax [2{x + t)+a] 



13/2 



V* = \ + Sax - f- [2(x + t) + af 2 + ^ [2(x - t) + af> + ^ 
Notice that (p cl coincides with (p b \ from formula (37). 
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Thus, the point x 3 (t) (see Fig. 8) accelerates from the initial velocity +2 to 
i 3 (3a) = 9/2. At the time t = 3a the points x ,xi meet each other, and the 
support of cp c i vanishes. It is the beginning of the fourth stage of evolution of (p. 



When considered in the context of the interacting strings mentioned in the 
Introduction, the solution discussed above provides an interesting picture of first 
stages of merging of the strings. In the initial state at t = 0, frame (a) in Fig. 9, 
the strings form a 'Y' shape with the vertex at the point x — 0, and the angle a 
between the two rectilinear pieces of strings in the region x > 0. In the region 
x < the strings have already merged. In the time interval [0, a) the vertex 
moves with the velocity -1 until the cusp is formed, see frame (c) in Fig. 9 - this 
happens at the time t = a. At this moment the bubble appears, which grows 
and moves along the strings in accordance with the self-similar solution of the 
signum-Gordon equation. Such smooth evolution goes on until t = 2a. In the 
third time interval, i.e. when t 6 [2a, 3a), the bubble growes and still moves 
along the strings, but now it is distorted by two light fronts xi(t), X2(t) travelling 
along it, as follows from the solution cp c given by formula (38). The evolution can 
be calculated also for times t > 3a, but we shall not dwell on it. 
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X A 




(a) (b) (c) (d) (e) 

Figure 9: The process of merging of two strings as described by solutions (34), 
(37). The 'bubble' visible in frames (d), (e) is a part of the self-similar solution. 
Frame (c) shows the cusp formed at the time T\ — a. Frames (d), (e) show snap- 
shots of the two strings at the times ti,t 2 such that a < t\ < t 2 < 2a. 

5 Summary and remarks 

1. The signum-Gordon equation appears in several problems of classical physics. 
It has the peculiar discontinuous nonlinearity, but this does not mean that the 
equation is intractable. We have constructed the full set of its self-similar solutions 
with initial data (5). The manifold of these solutions, presented as the map in Fig. 

1, is amazingly rich. Equally astonishing is the fact that it has been possible to 
find exact analytic forms of all solutions. 

The self-similar solutions may also appear as a component of generally non 
self-similar solution. This possibility is illustrated by the example presented in 
Section 4. It is interesting that the pattern of evolution of does not depend on 
the slope a in the initial data (33) - probably this property reflects the scaling sym- 
metry of the signum-Gordon equation. 

2. There are several obvious directions for extending our results. First, one may 
consider the four-parameter set of self-similar initial data (4) instead of the two- 
parameter set (5). It seems that not always the corresponding solutions can be 
obtained by straightforward combining solutions we have given above. Therefore, 
the map of such solutions will be truly four-dimensional. 

One may also be interested in stability of our solutions. For that matter, we 
have not seen any instability in the numerical investigations of the self-similar 
solutions. This seems to indicate that unstable modes, if present at all, grow rather 
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slowly. 

Probably the most interesting topic is related to dynamics of topological com- 
pact solitons (compactons) found in Q. These solitons have a piece- wise parabolic 
shape, and the pertinent field-theoretic model contains sectors which are described 
by the signum-Gordon equation. Therefore, one may expect that certain self- 
similar waves will appear in such processes as scattering of the solitons, or relax- 
ation of a single excited soliton. Perhaps one can even provide the exact analytic 
description of certain stages of such processes. This would give the much de- 
sired counterbalance to purely numerical investigations which are dominant in the 
literature so far. 
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